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Abstract 

We present in a local form the time dependent effective description of a su- 
perfluid Fermi liquid which includes Landau damping effects at T 7^ 0. This 
is achieved by the introduction of an additional variable, the quasiparticle 
distribution function, which obeys a simple kinetic equation. The transport 
equation is coupled with first order equations for the Goldstone mode and 
the particle density. We prove that a main feature of this formulation is its 
Hamiltonian structure relative to a certain Poisson bracket. We construct 
the Hamiltonian to quadratic order. 

1. Introduction 

The effective description at finite temperature of the broken symmetry 
phase of a Fermi gas in the BCS-BEC crossover is complicated by the ef- 
fect of Landau damping, which results in a highly nonlocal time-dependent 
Ginzburg-Landau theory [ij. Apart from the fact that the derivation of such 
an effective Lagrangian at T 7^ for either the complex order parameter or 
only the Goldstone mode is rather tricky even at quadratic level 
this nonlocal formulation makes it difficult to simulate numerically the real 
time evolution of non-equilibrium processes, such as oscillations of trapped 
Fermi gases. 

Indeed, a similar situation arises in dealing with collective effects and 
dynamical screening of Abelian and non-Abelian plasmas at high tempera- 
ture or high density. To leading order in the coupling constant, consistency 
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requires the inclusion of a set of one-loop diagrams termed "hard thermal 
loops", which are derived from a nonlocal effective Lagrangian Fortu- 
nately, the equations of motion can be written in a local form by introduc- 
ing auxiliary fields js], lof. In the Abelian case, the auxiliary field obeys a 
linearized Vlasov equation corresponding to the coUisionless regime in the 
plasma. Remarkably, it turns out that the resulting equations form a Hamil- 



tonian system with a noncanonical bracket structure [10|, lU 

In this paper, we turn to the question of the derivation of the low-energy 
dynamics of the phase of the order parameter in the superfluid phase, in 
connection with the above analogy. We present a simple derivation of the 
linearized equations of motion in a local form. Moreover, we prove that 
this local formulation is Hamiltonian, and H and the Poisson structure are 
completely identified. 

2. The equations of motion 

We begin by introducing the action for the system. In terms of the Nambu 
spinor = ijj^) the action for the two-component balanced Fermi system 
is written as 



2m 



S = j dx j dt^H^X) 

+ r+A(X) + r-A*{X)] ^(X) - J dx j dt y^*^, (1) 

where = | (t^ ± ^t^), and the r-^ 's are Pauli matrices. Here n is the chem- 
ical potencial, g\ is the bare coupling parameter, and we have included an 
arbitrary trapping potential V^xt(-^)- The complex field A(X) performs the 
Hubbard-Stratonovich decoupling of the quartic interaction between fermions 
in the BCS channel. 

In order to derive an effective theory when the symmetry [/(I) is broken to 
Z2, it is convenient to express the above Lagrangian in terms of a Goldstone 
field e{X) and "heavy" fields, ^(X) and A(X) to be integrated out 

^(X) = e'"'^(^)$(X), (2) 
A(X) = e2*^(^)A(X). (3) 

Let us choose the non-zero expectation value of the scalar field as real, (A) = 



2 



Aq. The condition that A(X) does not contain the Goldstone mod^ turns 
out to be Im ^A(X)*Aoj = 0, so the heavy field A(X) needs to be real with 
arbitrary sign. In terms of the covariant derivative 

D^^ix) = d^^ix) + td,,e{xy^{x), (4) 

the Lagrangian becomes 

2m \ J gA 

and the Noether current of the U{1) symmetry, which now is realized as 
= 6 A = and S9 = e, is given by 

= n = ^V^^, (6) 

jk = ^ f^t^^ _ 4$t$) + $V35 (7) 
2mi V / m 

To effectively integrate out the fermionic degrees of freedom we note that 
the total Hamitonian for these may be viewed as H = Hq + H^xt, where 

Ho = J dx^^ (-^'^ - - r'Ao^ ^, (8) 

and 

+ ^ f^tv$ - V$^$V V0 + 5"^^^^ I • (9) 

2mi\ J 2m J ^ ' 

Therefore i^ext couples the system to an applied perturbation given by ct{X) = 
A{X) — Aq, the gradients of 9, and Kxt- Following the procedure reviewed 
in Ref. [isj in the framework of the high temperature regime of QCD, the 
idea is to compute the induced changes in the expectation values of n, J 



^For a detailed discussion on how to construct effective Lagrangians in the case of 
spontaneously broken symmetries, see Weinberg [l^ . 
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and the pairing field \^'''r^\^ by the external perturbations {d^9,a,Vcxt}- 
By ignoring nonlinear corrections, we need the retarded response functions 
Xab{X) = —iii^i^): B{0)])9{t) whose Fourier transforms are denoted by 
XABiQY In the linear response approximation the induced changes take the 
form 

= -Xnn{Q)v{Q) + Xnj{Q)^q''0{Q)-Xnl{QMQ), (10) 

s{J'm = -x'}niQ)v{Q) + xyj{QWe{Q)-x'}i{QMQ), (n) 

5(M/ViM/) = -Xln{Q)v{Q) + Xu{Q)^q'0{Q)-Xn{QMQ), (12) 

where 

v{Q)=tooe{Q)-V,,t{Q). (13) 

From a functional point of view, if F*^^) [9, a] is the quadratic approximation 
for the effective action that one obtains after integrating out the fermions, and 
S[a] = —g]^^ J cr'^ is the bare action, then the integration of the heavy field a 
is simply gaussian. This produces an effective action for the Goldstone mode 
of the form F*^^) [9, a] + S[a] — | In det (— 2(yf^^ — xii) , where a is the solution 
to the saddle point condition 5(5'[cr] + T^'^^9, a])/Sa = 0. The last term, 
related to the determinant, does not depend on 9, but gives rise to important 
corrections to the thermodynamic properties evaluated in the mean field 



approximation [15|, Il6|. However, for the purpose of deriving an effective 
description only in terms of the phase 9, such a term will be ignored. Thus, 
in this approximation, the a field may be effectively integrated out by simply 
adjusting its value to the solution of the gap equation 

^ = 5($Vi$). (14) 

fl'A 

Exploiting the fact that n and 9 are canonically conjugated variables, and 
imposing conservation of the Noether current once a has been eliminated, we 
may derive a set of equations for the time derivatives of 6{n) and 9. These 
equations encode the effective dynamics. 

The explicit expressions of the response functions are easily computed 
at small frequency and momentum. By keeping the 0(A°) terms after the 
scalings co — )■ \co, q — )■ Xq, the general form of Xab{Q) when T 7^ is the 
sum a regular piece independent of Q, and a contribution due to Landau 
damping which is non- analytical at Q = 0: 

XAsiQ) = Xab + 2 / — — ^Fab — — ^r^- 15 
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The non- vanishing regular parts are given by 



Xni = = - y (^-^ tanh (17) 
xSt^ - -r/^4tanh^, (18) 



with the standard notation, E, = ^/^l + Aq, = k"^ /2m — /x. The form of 
Xjj is entirely due to the last term in Eq. ([7]). On the other hand, the factors 
in the integrand of the Landau damping contributions are given by 



Fnn{k) 


a 


Fni{k) 


= F^k) 


Fll{k) 


EV 


F'jjik) 


k^k^ 


FLik) 


= F^^j{k) 


F'nik) 


= F,!j{k) 



Ao6 
El 



mEk 
mEk 



(20) 
(21) 
(22) 

(23) 
(24) 
(25) 



These last results determine the nonlocal part of the departure from the 
equilibrium values of the quantities of interest. For instance, the nonlocal 
contribution to 6{n{X)) in space-time is obtained by combining the above 
terms according to Eq. ( fTOl) 

X/^^/n^\ nonlocal r, f Q ' kEk n'p{Ek) ,x , 
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where 

6EiQ,k) = ^-§^v{Q) + ^^a[Q) + rk.q'M 

= + (27) 

In order to present the equations in a local form, it is natural to introduce 
a new variable w{Q; k), defined as 

(28) 

which measures the departure of the distribution function from equilibrium 
when T 7^ 0. From this definition and Eq. fl26l) it follows that 

[u-q- VkEk) w{Q; k) = -n'p{Ek)VkEk ■ q 6E{Q, k), (29) 

and the equation of motion for distribution function w{X; k) takes the form 
of a transport equation without collision term resembling a linealized Vlasov 
equation 

{dt + VuEu ■ V,.)^(X, k) = n'p{E„)VkEu ■ VJE{X, k), (30) 

where SE is the induced change in the energy of the quasiparticle due to the 
applied perturbation 

SE{X, k) = -^v{X) + ^a{X) + k . (31) 

Thus, by combining the expressions for FAsik) with the definition of w, one 
finds the total changes in a local form 

__^^(X;fc), (32) 
5{J\X)) = + 2 [^^^w{X-k), (33) 

5(^t^i^(X)) = -xiVaiX) - Xniv{X) - 2 I ^,^^^{X- k). (34) 



The transport equation ( l30l) was derived long time ago by Betbeder- 
Matibet and Nozieres ^3\, and more recently by Urban and Schuck jisj. In 
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these approaches the starting point is the Heisenberg equation of motion for a 
matrix distribution function, whose diagonahzation leads to the above kinetic 



equation. As Leggett [19|, 120[ and Betbeder-Matibet and Nozieres [17[ have 
shown, it is possible to consider Fermi-liquid effects by adding to SE{X, k) 
an extra term, "^f,, fkk'w{X,k'), where fkk' describes interactions of two 
elementary excitations. 

To eliminate the a field we use the gap equation fll4l) . The relation 
between the s-wave scattering length and the bare coupling constant 

m I d?k I 

(35) 



4vras J {2n)^ 2ek 

together with the gap equation for Aq at T 7^ 0, 



m 



tanh^-1), (36) 



2™, J (27r)3 \Ek 2 efc 

produce = —2g]^^ — Xnn- Combining this result with Eqs. (fT4l) and ([3 
we obtaiiu 

Xni 2 r d^k A 



Finally, the conservation of the particle number 

dtS{n{X)) + V ■ 6{J{X)) = 0, (38) 

yields the equation of motion satisfied by v, thus completing the set of dy- 
namical equations for the variables {w{X; k), 9{X),v{X)}. This is found to 
be 

- ^^ni\dv f d^k ( ik Xni Ao\ _ „ _ ,^ 



XlJ dt J (27r)3 XnnEk, 

J (27r)3 \Ek XnnEj ^> ^ ^ ^\ m J 
+Mv2^ + 2 / tI^- ■ ^MX, k) = 0, (39) 



The term proportional to v was omitted in 17 1 
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where we have used Eqs. (15U]) and (157|) . as well as the vanishing of the angular 
integration J dflk^kEkfik) for any isotropic function f{k) such as ^k/Ek and 
Ao/Ek- Thus, the dynamical equations fl30l) . fl39|) . together with 

|^ = -t;(X)-V;xt(X), (40) 

and the equation fl37|) for cr(X) form a closed set of local equations for the 
effective low energy theory of the superfluid Fermi liquid at the one-loop 
level, when Fermi liquid effects are ignored. 



3. Hamiltonian formulation 

As we have pointed out before, a similar situation to the one posed by 
the Landau damping terms in the present context also appears when one 
attempts to derive a local time-dependent effective Lagrangian for the soft 
degrees of freedom of a gauge theory at high temperature or density. In that 
case, at the expense of introducing a new kind of degrees of freedom, one 
can reformulate the equations of motion as a system of local equations. The 
new variables represent the fluctuations of the charged (or coloured) particle 
distributions, and satisfy kinetic equations with external and induced fields 
due to these fluctuations It turns out that the complete set of dynamical 



equations are Hamiltonian with respect to certain Poisson brackets |10l . Ill 
It is natural to ask whether a similar formulation can be given in this case, 
and if so, what the Hamiltonian structure describing the low energy theory 
of the Fermi superfluid would be. 

To address this question, it is convenient to replace the variable v{X) by 
the more natural ni{X) = 6{n{X)), and to consider {w{X; k), 9{X), ni{X)} 
as the set of dynamical variables. This choice exploits the fact that the par- 
ticle density and the Goldstone mode are canonically conjugated, as follows 
from the role of the particle number operator as the generator of the U{1) 
symmetry, and the inhomogeneous transformation law for 6, 66 oc i[N,6]. 
Therefore the Poisson bracket of these variables may be written as 

{m{t,x),e{t,y)} = 6{x-y). (41) 

The remainder Poisson structure may be guessed by noting that a kinetic 
equation of Vlasov type usually takes the Hamiltonian form 

dJ{t,x,k) = {f,H}. 



8 



Here the relevant bracket is the Poisson-Vlasov bracket 



21 



{m,G[f]} = I ^'^^'^/{^'^} ^' (42) 

where {, is the canonical bracket on the single particle phase space 
spanned by {x,k), 

{/, g}^^ ^ VJ ■ Vkg - V^g ■ V,/. (43) 

The functions f,g are distribution functions on M^. Since we are interested 
in the linearized equations around equilibrium, where / = np^E^) +w{x, k), 
the Poisson structure will be chosen as 

{F,G} = pjd'xJd'knpiE,)!^^,^^^ 
/•^3 f5F5G 6G6F\ 



The first term corresponds to a frozen Lie-Poisson bracket [2l|] which is ob- 
tained from (l42l) by linearization; the value of the factor P will be determined 
shortly. The Jacobi identity for the frozen bracket follows from the general 
properties of Lie-Poisson brackets [2T]. The second line corresponds to the 
canonical Poisson structure of Eq. fHTl) . 

It remains to see that the equations of motion we have found before are 
Hamilton's equations 

dtW = {w,H} = Pn'p{Ek)VkEk-V,(^j-^^y (45) 

- 

= {n,,H} = ^^^, (47) 
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for some quadratic functional H[w^9^n^. With the aid of the relations 



'X.nn ~^ Xnl 



+ 2 



d^k ( Xni Ac 



(2vr) 



a(X) 



Xr. 



-f'fc Xnn Ek 
Xnl 



w{X,k) 



(48) 



Xr 



-ni(X) 



+ 2 



Ao _^ Xnl ^fc 



(27r)3 XnnEk 



w{X,k) 



(49) 



we may express the change in the energy of the quasiparticle in terms of the 
new variables 



6E{X, k) 



Xr. 



6 Xnl Ac 



Xnn ~l~ Xnl \Ek Xnn E]^ 



ni{X) 



+k.^^^ + 5E('\X,k) 
m 

where SE'^^'^ denotes the contribution linear in w 
6E^'\X,k) -- 



(50) 



2Xnn 


6 


r <fk' 




Xnn ~^ Xnl 


EkJ 


(27r)3 


\Ek' 



Xnl Ac 



Ao f d?k' 



Ao _^ Xnl ik> 



w(X, k') 



w{X, k') 



Ek J (27r)3 yEk' ' XnnEk' J "^^^^'-^J • 

By comparing Eq. P7|) for (9tni and Eq. f l33|) , it follows that the Hamiltonian 
must contain exactly the term 



On the other hand, its derivative with respect to w when inserted into the 
RHS of Eq. psj) produces 

W(X)\ 



Pj^n'p{Ek)VkEk ■ ■ 



so, in order to match the corresponding piece of dfW in Eq. ( 130|) . P must 
assume the value P = (27r)^/2. Now, since 
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6w{X, k) 



d'y d\ 6E^'^ (t, y, q)w{t, y,q) = 2 6E^'^ (X, fc) , (52) 
10 



one can check immediately by functional derivation that the Hamiltonian 

Xnn ni{Xf 



H[w,e,ni] = j (fx 



Xnn ~^ Xnl 



2m 

3 f d^k 2Xnn f Xnl ^0 



+ M'M 3 -2 T-2 f- - n,{X)w{X, k) 

{27t)^ Xln + Xnl \Ek XnnEk 



+ / d^x 



.^44^ + '^^^'H^, k)w{X, k)] ,(53) 



d^k 

and the bracket flH|) yield the equation of motion for {w, 6,ni}. Although the 
Poisson structures are decoupled, the effective Hamiltonian contains terms 
mixing all variables. As n'p{Ek) is a monotonic function the regularity of the 
integrand is guaranteed. 

When Aq — )■ the coefficients Xnn and Xni vanish, and the most singular 
terms of the Hamiltonian are grouped as 

X ^ [ f f d^k \^ 

To keep the energy finite, the change in the particle density for vanishing 
x's is restricted to ni = 2 J^sign(^fc)w, according with Eq. (132|) . Thus the 
equation of motion for 9 yields dtO = 0, and this variable in no longer time 
dependent. Now the Hamiltonian becomes 

d^k w{X,kf 



(27r)3 n'p{ik) 
d^k 



+2 J d-'xVe.tiX) J -^sign{^k)w{X,k) 

+^/^^^/(0^-^-^^'^)' (^^) 

where we have used n'p{\^k\) = n'p{^k)- The equation of motion for dtUi 
obtained by variation of 6 is not an independent equation as it corresponds 
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exactly to the integration of 2 sign{^k)dtw. If the initial condition is cho- 
sen 6{t = 0,x) = constant, the V^(a?)-terms in the Hamiltonian may be 
ignored, and one recovers the Hamiltonian of a non-interacting Fermi gas in 
an external potential [§. By contrast, in the limit of zero temperature when 
n'p{Ek) — > 0, the distribution function w must vanish in order to keep the 
energy finite. Therefore one recovers a description in terms of rii and 9 alone. 

4. Conclusion 

To conclude, we have provided the Poisson structure and an effective 
Hamiltonian for the low energy description of a superfluid Fermi gas in the 
collisionless regime. The basic dynamical variables are the quasiparticle dis- 
tribution function, the Goldstone mode, and the particle number density. 
Apart from the limitations inherent to the linear approximation to the equa- 
tions of motion we have made, the most serious limitation of this approach is 
that it neglects Fermi liquid effects. A more accurate treatment would require 
the addition of the contribution 1/2 j dk dk' fk.k'w{X, k)w{X, k') in terms of 
the appropriate Landau parameters to the above effective Hamiltonian. 
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